Hybrid configuration content of heavy S-wave mesons 
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We use the non-relativistic expansion of QCD (NRQCD) on the lattice to study the lowest hybrid 
configuration contribution to the ground state of heavy S-wave mesons. Using lowest-order lattice 
NRQCD to create the heavy-quark propagators, we form a basis of "unperturbed" S-wave and 
hybrid states. We then apply the lowest-order coupling of the quark spin and chromomagnetic field 
at an intermediate time slice to create "mixed" correlators between the S-wave and hybrid states. 
From the resulting amplitudes, we extract the off-diagonal element of our two-state Hamiltonian. 
Diagonalizing this Hamiltonian gives us the admixture of hybrid configuration within the meson 
ground state. The present effort represents a continuation of previous work: the analysis has 
been extended to include lattices of varying spacings, source operators having better overlap with 
the ground states, and the pseudoscalar (along with the vector) channel. Results are presented 
for bottomonium (T, n b ) using three different sets of quenched lattices. We also show results for 
charmonium (J/ip, r] c ) from one lattice set, although we note that the non-relativistic approximation 
is not expected to be very good in this case. 

PACS numbers: 11.15.Ha,12.38.Gc 



I. INTRODUCTION 

The existence of valence gluons in bound quark sys- 
tems, a theoretical possibility considered for a long time 
in QCD, continues to elude confirmation. Allowance 
for gluonic excitations increases the range of possible 
hadronic quantum numbers (J PC ) beyond those pre- 
dicted by constituent quark models. Exotic states should 
appear (for recent reviews see Refs. 0,0) an d, in fact, 
one such state (1~ + ) has been observed recently 0, but 
the underlying structure of the state (hybrid meson, four- 
quark state, meson molecule, etc.) has not been deter- 
mined. When considering gluonic constituents, however, 
the resulting state need not be exotic; the valence glu- 
ons may also combine with the quarks and antiquarks to 
form a state which may otherwise be formed without the 
gluonic presence. In this non-exotic scenario, a meson 
(or baryon) state should consist of a mixture of configu- 
rations; not only the case where only the valence quarks 
and antiquarks appear, but also those where a gluonic 
excitation is present: a hybrid configuration. 

As a relevant example - since this is one of the systems 
we study in the present work - we may look at the vector 
meson state for bottomonium (T, 1 ). We may envi- 
sion the ground state for this system as a bottom quark 
and antiquark in a color singlet, a relative S-wave, and a 
spin triplet. However, the true ground state should also 
have a contribution where the quark and antiquark are 
in a spin singlet and a color octet, the spin of the meson 
and the overall color singlet being ensured by the gluonic 
excitation: 

|T) = A s \bb) + A h \bbg) 

= cosO\lS{l~~)) + wa.9\lH{\~~ )). (1) 

It is just this type of (albeit simplified) two-state system 



we consider in the present work. We also consider the 
h heavy S-wave meson: 



|%) = cos0 / |lS((n + )) +sin0'|Lff((T + )). 



(2) 



Hybrid-quarkonium configuration mixing has been 
considered before in the framework of the MIT bag model 
4 , and with the use of an adiabatic potential model 
6 1 . We compare results from this on-going lattice work 
LULU with these previous results. 



II. LATTICE METHOD 

We work in the heavy-quark limit, so we use the 
non-relativistic expansion of lattice QCD fit! TTTL IT^ . 
To evolve our quark propagators we use a time-step- 
symmetric form for the transfer matrix 0, Il3| : 



(x, t + a) 
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t+a \ / t 

x(l-S t ', t a5H) t (f>(x,t). (3) 



The Hamiltonian applied at all time slices, 7io, accounts 
for the kinetic energy of the heavy quarks: 



Ho = 



A 2 
2m„ 



(4) 



where the A 2 is the lattice covariant Laplacian. At 
one intermediate time slice, t', the lowest-order spin- 
dependent interaction, 



SH 



' 2m, 



-a-B, 



(5) 
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is applied, thereby allowing a spin flip of the quark (or 
antiquark) in exchange for the emission or absorption of 
a gluonic excitation; i.e., configuration mixing. 

The local value of the chromomagnetic field is calcu- 
lated using the clover formulation, averaging the fields 
generated from the four plaquettes surrounding a lattice 
site: 

%fa) = ~[Ufa)Ufa + fi)Ufa + i>)Ufa) 

+ U v {x)Ufa + 0- fi)Ufa - (£)Ufa - A) 

+ Ufa - fi)Ufa — i> — ji)Ufa — v — jj,)U u {x - v) 

+ Ufa - v)Ufa - v)U v (x -v + (i)Ufa)] , (6) 



^fa) = ~ [^u(x) - n^fa)] - ilmpifVOs)]- (7) 

The chromomagnetic field arises from the spatial compo- 
nents, 



Fjfa) = s jk igBi(x) 



(8) 



Tadpole improvement [T^ is also included, the factor uq 
being calculated via the average plaquette and applied to 
all the link variables: 



Ufa) 



Ufa) 

Uq 



3~ \iv\X) 



(9) 



Creating the heavy-quark propagators with this form 
for the evolution operator, we then use appropriate op- 
erators at the source and sink time slices to project out 
the desired meson states. The heavy-meson operators 
we use, along with the corresponding quantum numbers, 
may be found in Table 1 of Ref . . 

Two different types of quark sources are used: a ran- 
dom wall (RW) and a Coulomb-gauge-fixed wall (CW). 
The RW source is an incoherent collection of point 
sources, the average meson propagator having contribu- 
tions only from where the quark and antiquark start at 
the same location. The CW source provides a "maxi- 
mally smeared" source, with contributions from spatially 
separated quarks and antiquarks. In both cases, the sink 
end is simply a sum over points where both the quark 
and antiquark coexist. 

We fit the S-wave propagator with the following form: 



C s (t) = A ls e- m ^ + A 2s e 



-m 2s t 



(10) 



allowing a determination of the 2S-1S mass splitting. We 
also have P-wave correlators, which we fit (at least for the 
CW source) with only a single mass, 



C p (t) = A 2p e~ 



(11) 



We use the 2P-1S mass splitting to set the lattice scale. 
A single-mass fit is also used for the hybrid correlators, 



C h (t) = A lh e 



(12) 



All fits to the propagators use the full covariance matrix 
to account for correlations among the different Euclidean 
times. 

After application of the spin-dependent term in the 
Hamiltonian at t' (or t"), a signal appears for the 
"mixed" correlator, hybrid — > S-wave (or vice versa). We 
fit these correlators in the region t > t' (or t > t") with 
the forms: 



and 



C^(t',t)=A lhs (t')e- m ^ t - t '^ (13) 



C^(t",t)=A lsh (t")e- m ^ t - t "\ (14) 



Looking more closely at the amplitude for the first corre- 
lator (hybrid — * S-wave) , we can reason that there should 
be factors from the overlap of the source operator with 
the hybrid, the overlap of the sink operator with the S- 
wave, the exponential decay of the hybrid state before t' , 
and the matrix element with which we are interested: 



41/2,1/2 



lS\aSH\ lH)e 



-m lh t 



(15) 



Knowing the masses and amplitudes from the standard S- 
wave and hybrid correlators, we solve for the off-diagonal 
matrix element of our two-state Hamiltonian. This is re- 
peated for larger values of t' to find a plateau in the final 
result, where we may be sure that only the ground-state 
contribution from the source appears. A similar proce- 
dure may be followed for the S-wave — > hybrid correlator. 

There is a complication, however, which arises for our 
CW-source correlators: the operators at the source and 
sink ends are not the same. At one end we have a 
Coulomb-gauge-fixed wall (cw) source, while at the other 
end there is a point (p) sink. The amplitude for the S- 
wave correlator therefore has the form 



while that for the hybrid is 



,1/2 .1/2 



(16) 



(17) 



From each of these products, the mixed correlators in- 
clude only one of the factors (rather than just the square 
root of each amplitude) : 



1 1/2 A 1/2 /KU^nm ,,-mut' 



A lhs (t') = A^ p A^ cw {1S\oSH\1B), 



and 



A ls fa')=A{ / h 2 A 1 1 / s l(lH\aSH\lS), 



-mi s t" 



(19) 



In order to get the appropriate cancellations of ampli- 
tudes, we thus need to use a geometric mean: 



|(15|o(SW|lJ?)| 



1 A lhs {t')A lsh {t") 
Ai s A lh 



(20) 



TABLE I: Quenched NRQCD runs. 



p 


x N t 


wo 


am q 


n 


sources 


# configs. 


7.75 


16 3 x 32 


0.8800 


3.2, 3.6 


2 


RW,CW 


220 


8.00 


20 3 x 64 


0.8879 


2.5, 2.8 


2 


RW,CW 


240 






n 


0.7, 0.8 


3 


CW 


170 


8.40 


28 3 x 96 


0.89741 


1.8, 2.0 


2 


RW,CW 


76 



at large t' , t" . 

We average the correlators over sets of quenched lat- 
tices, generated using a Symanzik 1-loop improved gauge 
action [3 [lj| E [TJ Eg- The relevant parameters for 
these runs are listed in Table [I] 

To set the physical scale for the lattice spacing and 
determine the physical quark mass, we use the same pro- 
cedure as was used previously for bottomonium spec- 
troscopy with NRQCD 19]. For the lattice spacings, 
we use bottomonium mass splittings, such as the spin- 
averaged 2P-1S mass difference (A% - My = 440 MeV) 
and the 2S-1S difference for the T (My - My = 563 
MeV). We also create non- zero-momentum S-wave cor- 
relators and use the resulting dispersion relations to de- 
termine the kinetic masses of the ground states. Using 
the results from two quark masses, an interpolation (or 
extrapolation) may be made to match the kinetic S-wave 
mass to that of the T {My = 9.46 GeV), thus arriving 
at a physical bottom quark mass. 

In an attempt to determine the radiative correction c B , 
we also create correlators with the spin-dependent term 
applied at all intermediate time slices and with the val- 
ues of c B = 1 and 2. We then determine the resulting 
S-wave hyperfine splittings, which to lowest order should 
be quadratic in this term: mi s (l ) — rni s (0 '") oc c B . 
There is, however, a danger in assuming this to be 
the only (or the largest) contribution to the hyperfine 
splitting: previous investigations of quarkonium spin- 
dependent splittings in lattice NRQCD [2(], 0] display 
significant contributions from various 0(i> 4 ) and 0(i> 6 ) 
terms in the velocity expansion (and poor convergence 
of this expansion for cc systems). Therefore, a value of 
c B (a) determined in this way includes not only the de- 
sired radiative correction, but also systematic effects due 
to the neglect of other relativistic terms in the NRQCD 
expansion. To further complicate the matter, the 
bottomonium state (rj b ) has yet to be observed experi- 
mentally. 



III. RESULTS 

The Coulomb-gauge-fixed wall (CW) sources provide 
good overlap with the desired meson states and rea- 
sonable fits are obtained earlier than for the random 
wall (RW) sources. This can be easily seen in effective 
mass plots of the P-wave (Fig. ^) and hybrid correla- 
tors (Fig. |2J) • The CW sources provide correlators which 
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FIG. 1: Effective masses for the S-wave (1 ) and P-wave 
(0 ++ ) correlators from the (5 = 8.0 lattices with am, = 2.5. 
The dashed horizontal lines denote the masses (am lg , am 2p ) 
from the (correlated) fits to the unperturbed CW-source cor- 
relators, along with the corresponding fit ranges. 



approach plateaus earlier than their RW counterparts. 
The CW correlators are less noisy as well; this is espe- 
cially noticeable for the hybrids. Both S-wave correlators 
take a substantially longer time to approach a consis- 
tent plateau. These are reliably fit by the two-mass form 
found in Eq. (|10fl . Also included in the figures are exam- 
ples of the mixed correlators. In Fig.^the effective mass 
for the H— >S correlator is shown for the region t > t' = 5. 
Although the mixed correlator is much more noisy, the 
agreement with the mass from the unperturbed S-wave 
is good and we may fit this correlator over a large range 
of t to extract the mixed amplitude, Eq. I)18[l. In Fig. [3 
we also show the effective mass for the S— >H correlator 
for t > t' = 12. Here we see the difficulty which arises, 
namely the quality of signal, for the S— >H correlators in- 
volving the point-like sinks. Nevertheless, for each value 
of t' , we have at least a few significant points which give 
masses consistent with those from the unperturbed CW- 
source hybrid correlator. We use these to extract the 
corresponding mixed correlator amplitudes, Eq. l|19fl in 
this case. A complete list of the fits used for our results, 
along with the amplitudes, masses, and \ 2 values may 
be found in Ref. Q- 

As an aside, we point out the significant curvature in 
these plots for the H— >S and S— »H correlators in the 
region t' w t, even though the excited-state contribu- 
tions from the sources should be minimal: e~ m2ht and 
e ~m 2s t ^ respectively, become small. These are most 
probably mixings with higher-mass configurations: e.g., 
(25|a<57<| IH) and (2H\a5H\ IS). In principle, the ex- 
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FIG. 2: Effective masses for the _+ hybrid correlators from 
the j3 — 8.0 lattices with am q — 2.5. The dashed horizon- 
tal lines denote the mass (am lh ) from the (correlated) fit to 
the unperturbed CW-source correlators, along with the cor- 
responding fit ranges. 



TABLE II: Lattice spacing determinations with c B — 0, CW 



0.344 fm 



-1/2 



440 MeV 



source. Physical scales used: n 
■ M %b - M r = 440 MeV ; My - Mr = 563 MeV. The 
denote values used for a throughout this work. 
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— amq — 


physical scale 


n^ 1 (MpV\ 

Ui t Ivlc V j 


7 7^ 
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Tl/Q — 








DC 


a 2 n = 


D 1 6^9/471 


1 089('1 ^ 




3 9 


a \ M Xb 


— IVlr ) — \).o£o\Lo) 






6 . 


n( M- 

a \ lvl Xb 




looo^ol ) 




'< 9 


ai^Mr' 


— /If-vl — 487f'60 s l 






3.6 


a( Mr> 


— A/f-v"! — 434('34~) 


i soon oo^i 


s nn 


X. 


Ti/a = 








DC 


a 2 K = 


OQQ55H 0^ 






9 r i 

Z .U 


a \ 1V1 Xb 








9 S 
Z.o 


a(M Xb 


- Mr) = 0.2431(40) 


1011 /or^ 




2.5 


a(Mr' 


- Mr) = 0.3233(88) 


1742(48) 




2.8 


a(Mr> 


- Mr) = 0.3131(73) 


1797(43) 


8.40 


CO 


n/a = 


3.7301(69) 


2136(4) 




CO 


an — 


0.04989(46) 


1970(9) 




1.8 


a{M %h 


- Mr) = 0.1749(50) 


2516(71)* 




2.0 


a(M Xb 


- Mr) = 0.1724(47) 


2552(69) 




1.8 


a(M r > 


- Mr) = 0.235(16) 


2400(160) 




2.0 


a(M r > 


- M T ) = 0.229(14) 


2460(140) 



traction of these configuration-mixing amplitudes should 
also be possible, especially with finer lattice spacing in 
the time direction (e.g., with anisotropic lattices p^ . l23| ). 
For our purposes here, however, we focus only upon the 
mixing between the lowest-lying configurations, IS and 
1H. 

The results for the lattice spacings are presented in Ta- 
ble [n] along with other determinations from the static- 
quark p otential using a modified Sommer parameter n 
|24ll25l ] and the string tension k. Looking at the NRQCD 
results, we are encouraged by the consistency (within 
the errors of at least one) of the 2P-1S and 2S-1S mass 
splittings. (We would like to point out, however, that 
our heavy-quark Hamiltonian leaves out many terms in 
the relativistic expansion and, for this reason, we do not 
claim to be presenting a very precise determination of 
the bottomonium spectrum; this has been studied by 
others using more elaborate forms of lattice NRQCD 
El EES-) There arc marked differences between the 
lattice spacings determined via the bottomonium mass 
splittings and those from the static quark potentials; in 
fact, the two static-quark determinations disagree. The 
string-tension results consistently give the largest lattice 
spacings (a in fm), while the bb mass splittings give the 
smallest. The main reason for the difference of these two 
extremes appears to be the quenched approximation. A 
comparison of the static-quar k p otentials for quenched 
and dynamical configurations |25[ has shown that, in the 
quenched case, the potentials do not display sufficient 
curvature: the Coulomb-like potential well at short dis- 



tances is not as deep as in the dynamical case and the 
linear string-like region is steeper. This should lead to an 
underestimation of the bb mass splittings in lattice units 
since these systems are relatively small. This would seem 
to explain the low values for the lattice spacings seen in 
Table [n] The string-tension results are thus expected 
to give the larger lattice spacings for the quenched lat- 
tices. (A cursory study |9( using the Born-Oppenheimer 
approximation and the quenched versus unquenched po- 
tentials [2{| supports these claims. Others [22L |26| have 
also found, with better energy resolution, discrepancies 
between 2S-1S and 2P-1S mass splittings on quenched 
lattices, suggesting the same effect.) In the spirit of pre- 
senting a self-consistent work, however, we use the lattice 
spacings provided by the 2P-1S bottomonium mass split- 
tings [23. 

In Table ITTT1 we present our results for the IS kinetic 
masses determined from the dispersion relations: 

^{V) = ^ w +mr, (21) 

Each case requires an extrapolation in quark mass to 
reach the physical value (Mff 1 = M T = 9.46 GeV). 

As described in the previous section, we use the ampli- 
tudes from the mixed correlators, along with the masses 
and amplitudes from the "unmixed" ones, to determine 
the off-diagonal element of our two-state Hamiltonian. 
Figures I3HH1 show the results for these matrix elements 
as functions of the time, t', at which the spin-dependent 



5 



0.15 



TABLE III: Kinetic masses for the IS states and the resulting 
(lattice-regularized) physical quark masses. The lattice scale 
is set using the spin-averaged 2P-1S mass differences found 
for the 66 systems. 



p 


am q 


aMff 1 


Mf! n (GeV) 


am b 


m b (GeV) 


7.75 


3.2 
3.6 


7.38(54) 
8.44(59) 


9.90(72) 
11.46(80) 


3.09(21) 


4.17(28) 


8.0 


2.5 
2.8 


5.509(77) 
6.178(90) 


9.87(14) 
11.19(16) 


2.41(4) 


4.34(7) 


8.4 


1.8 
2.0 


4.00(16) 
4.42(17) 


10.06(40) 
11.28(43) 


1.70(7) 


4.31(18) 










am c 


m c (GeV) 


8.0 


0.7 
0.8 


1.752(23) 
1.964(26) 


3.154(41) 
3.535(47) 


0.677(12) 


1.219(22) 



term is applied to the heavy-quark propagators. The 
CW-source, H^S results each have the remaining factor 
of 



.1/2 ,1/2 



1/4 , 1/4 



lh c 



A 1/2 A 1/2 



1/4 ,1/4 



A 

~A 



1/2 



(22) 



(the reciprocal for S— >H) which requires the use of a ge- 
ometric mean (see Eq. I2(J|1 to ensure its removal. The 
results chosen for the final geometric means are displayed 
with dotted symbols. All errors result from a single- 
elimination jackknifc routine. 

Figure [3] displays results for the configuration-mixing 
matrix element from one set of lattices (J5 — 8.0) and one 
value of the quark mass (am q = 2.5). Results are shown 
using both the CW and RW sources. The hybrid corre- 
lators for the latter, however, have their masses fixed to 
the values found from the CW-source hybrid correlators 
since these provide more reliable mass plateaus. The two 
horizontal lines show the ±ler limits for the geometric 
mean of the CW-source results. The hybrid^S-wave re- 
sults do not display very convincing plateaus, especially 
for the 1 channel. However, within the errors, the fi- 
nal results are consistent with the plateaus from the RW 
source. These results appear to be about 30% lower than 
previous determinations using only the RW sources 
thereby stressing the need for reliable hybrid mass deter- 
minations, which the CW-source correlators more readily 
provide (see Fig- EJ- 

Figures 0| and [5| display similar plots for the (3 — 7.75 
and 8.4 lattices, respectively. Here we see more convinc- 
ing plateaus for the matrix element from the hybrid— >S- 
wave correlators. 

In Fig. we show the results for one of the lighter 
masses (around m c ) for the f3 — 8.0 lattices. A quick 
look at the result for the 1 channel, however, shows 
that the non-relativistic approximation (at least at the 
level of simplicity in our heavy-quark Hamiltonian) may 
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FIG. 3: Comparison of the results obtained for the mix- 
ing matrix element via the CW source and the RW source 
with the hybrid mass fixed to that found with the CW hy- 
brid source. The dotted symbols mark the ones used in the 
geometric mean to get the CW-source results, the la ranges 
of which are denoted by the horizontal lines. 
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FIG. 4: The mixing matrix element vs the time of application 
of the interaction term. CW source, /3 = 7.75, am q = 3.2. 
The dotted symbols mark the ones used in the geometric mean 
to get the CW-source results. 



not be such a great idea: 

(IH\6H\1S) _ 0T4 
2m„ LT 



0.1 



(23) 



The mass required at this lattice spacing to reach 
the charm quark is also problematic for the expansion 
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FIG. 5: The mixing matrix element vs the time of application 
of the interaction term. CW source, /3 = 8.4, am, = 1.8. The 
dotted symbols mark the ones used in the geometric mean to 
get the CW-source results. 
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FIG. 6: The mixing matrix element vs the time of application 
of the interaction term. CW source, f3 = 8.0, am q = 0.7. The 
dotted symbols mark the ones used in the geometric mean to 
get the CW-source results. 



(am c < 1). In spite of these problems we carry on 
and present our results for charmonium, encouraging the 
reader not to forget the large systematic effects we in- 
troduce by neglecting other terms in our heavy-quark 
Hamiltonian and by simulating at such a small quark 
mass. 

All of the CW-source results for the configuration- 



TABLE IV: Results for the 1 S-wave/hybrid configuration 



mixing 


(CW 


source) . 






a 
P 


CLTTlq 


1 /I W InA'L/l 1 <?\ 1 
(ZD n\ ID / 1 




\\ Ln 1 1 /I 


7.75 


3.2 


0.0624(61) 


0.0566(27) 


0.058(3)c B (/3) 




3.6 


0.0582(61) 


0.0518(27) 




8.0 


2.5 


0.0566(33) 


0.0618(28) 


0.063(3)c B (/3) 




2.8 


0.0522(31) 


0.0564(27) 




8.4 


1.8 


0.0369(44) 


0.0618(36) 


0.064(4)c B (/3) 




2.0 


0.0352(44) 


0.0578(36) 




\(1H\J/1>)\ 


8.0 


0.7 


0.1349(84) 


0.1483(81) 


0.150(8)c s (/3) 




0.8 


0.1259(74) 


0.1393(70) 




TABLE V: Results for the 0" 4 


S-wave/hybrid configuration 


mixing 


(CW 


source) . 









am q 


\{lH\aSH\ 15>| 


[sin(0)l 




7.75 


3.2 


0.1032(58) 


0.0945(26) 


0.097(3)c s (/3) 




3.6 


0.0954(57) 


0.0862(26) 




8.0 


2.5 


0.0966(34) 


0.1031(29) 


0.106(3)c s (/3) 




2.8 


0.0887(32) 


0.0939(27) 




8.4 


1.8 


0.0767(60) 


0.1156(52) 


0.120(5)c B (/3) 




2.0 


0.0723(61) 


0.1074(53) 




\{lH\Vc)\ 


8.0 


0.7 


0.2327(77) 


0.2403(62) 


0.243(6)c B ()9) 




0.8 


0.2172(69) 


0.2270(56) 





mixing matrix elements, along with the corresponding 
mixing angles, are shown in Tables ITVl and Ivl The last 
column displays the final results after an extrapolation 
(linear) in quark mass to the bottom (or charm) mass 
determined previously (see Table II I II) . The hybrid con- 
figuration content of the ^ ground state is enhanced 
by the expected factor of V3 (due to spin statistics for 
single-gluon emission/ absorption 0,H^|) relative to that 
in the 1 channel. Up to this point we have ignored 
the radiative correction (c B = 1) and we point out the 
need for this factor, c B (a{0)), in the final column. To 
be precise, this factor should be included in the ma- 
trix element, but for small mixing angles, this correc- 
tion appears as a multiplicative factor in the final result: 
sin (9 sa (lH\SH\lS)/(m lh - m ls ). 

Meson correlators with the a ■ B term applied at all in- 
termediate time slices were also created in order to study 
the S-wave hyperfine splittings and to attempt to de- 
termine c B ((3). For these, the values of c B = 1 and 
2 were intended; however, a sign error was discovered 
(much too late) in our field-strength routine which forces 
us to work with c B = —1 and —2. (This has no effect 
upon our configuration-mixing results up to this point 
since they result from a single application of the spin- 
dependent term.) The resulting hyperfine splittings for 
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TABLE VI: IS hyperfine splittings. 



p 


am q 


„2 
C B 


71 /f 71 /T /1\/T~\A 

Mt — (MeV) 




c b(P) 


7.75 


3.2 


1 


16.7(1.3) 


(Mt 


- M % ) eX p/(17MeV) 






4 


63.6(5.1) 






8.0 


2.5 


1 


21.90(55) 


(Mt 


-M I) Jex P /(22MeV) 






4 


81.6(2.6) 






8.4 


1.8 


1 


25.2(1.1) 


(Mt 


-M % )ex P /(25MeV) 






4 


93.5(4.9) 






M J/4 , - M Va (MeV) 


8.0 


0.7 


1 


68.0(1.8) 




w 1.2 - 1.7 






4 


213.6(7.1) 







the S-wave states may be found in Table IVll 

Whereas our bottomonium systems seem to follow the 
c B proportionality and our charmonium systems do not, 
there is no experimentally measured mass for the h 
bb meson (jj h ) while the corresponding cc meson (r) c ) has 
been observed. So we report our findings for c 2 B (a((3)) 
from the bottomonium systems in terms of the Mt — M Vb 
mass splitting. For the charmonium results, we report 
two values for c B (/3). The larger value is the result of 
a linear extrapolation from through the c% = 1 point 
to M J/4 , - M Vc = 117 MeV at c| = c|(/3). If we be- 
lieve the physical quark mass shifts with the inclusion 
of the spin-dependent term and that this is the domi- 
nant effect in explaining why our cc hyperfine splitting 
is not oc c 2 B , then this is the result we would trust more. 
If, however, we believe there is a c B correction present 
and that this is the dominant effect, we should trust the 
lower value, which results from using an interpolation of 
the form a2C B -\-a^c B through the c B = —1 and —2 points 
and extrapolating to 117 MeV in the c B > region. If 
we use these results with that from bb at (3 — 8.0, we find 
M T -M Vb « 26-38 MeV and c B (f3 = 7.75) » 1.24-1.50, 
c B (j3 = 8.0) « 1.09 - 1.31, c B ((3 = 8.4) w 1.02 - 1.23. 
Any such determination of the radiative corrections, how- 
ever, is further complicated by the fact that it incor- 
porates systematic effects due to the termination of the 
NRQCD expansion [U|2l|. 

Given our previous worries about our cc results (u 4 ~ 
0.1; am c < 1), we may also consider the possible range of 
c B values using potential model values for the bb S-wave 
hyperfine splittings [3(1 : ^ 30 — 60 MeV. Assuming 
the a-B term to dominate this splitting results in c B {(3 = 
7.75) - 1.33 - 1.88, c B (/3 = 8.0) - 1.17 - 1.65, c B (fi = 
8.4) - 1.10 - 1.55; or roughly 

|(Lff|T}| ~ 0.076 -0.11 ; | (1-81%) | ~ 0.13-0.19 
and 

\(\H\J/%I))\ ~ 0.18-0.25 ; |(lfl"|r? c )| ~ 0.29 - 0.4. 

Each of these ranges lies between two corresponding val- 
ues (with and without the color Coulomb interaction) 



0.2 



55 0.1 







^^^^^^^ c 


|— -t— : 




$ 


O : 0" + 
❖ : 1"" 



a(M 2P - M 1S ) 

FIG. 7: Hybrid/S-wave configuration mixing angle vs lat- 
tice spacing. For each channel, the three lowest points (solid 
fit line) represent the (tadpole-improved) tree-level results: 
c B — 1. The middle three points result from setting the ra- 
diative corrections with the S-wave hyperfine splitting held 
constant at 30 MeV. The uppermost points result from the 
same splitting set to 60 MeV. The 1 points are shifted 
slightly to the left for clarity. 

determined within the MIT bag model Q. The upper 
limit of the J ftp result is consistent with that from an 
adiabatic potential model @. 

Again, however, the above determinations of c B (p) in- 
clude systematic effects due to the neglect of other terms 
in the velocity expansion which also contribute to the hy- 
perfine splittings [20I l2ll | . We are therefore left without 
precise values for the radiative corrections and may only 
state that if the a ■ B term does in fact dominate the 
S-wave hyperfine splitting, then these corrections appear 
to enhance the configuration mixing further. 

In Fig. \7\ we plot our results for the bottomonium 
configuration-mixing angle as a function of the lattice 
spacing. The plot includes the tree-level results, along 
with two others altered using radiative corrections given 
by assumed values for the S-wave hyperfine splitting due 
to the a-B term. Some positive features may be noted in 
this plot. First, the choice of the "physical" value for the 
hyperfine splitting does not seem to have a large effect 
upon the scaling displayed by the points (however, the 
values at any one lattice spacing are greatly affected). 
Also, at least for the ^ channel, the fixing of the hy- 
perfine splittings appears to slightly improve the scaling; 
the tree- level and radiatively-corrected values for the 1~~ 
appear to scale equally well. 

For the bottomonium systems on our lattices, where 
we have a > 1 / discretization effects are expected to 
be about the same order as those due to the relativistic 
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corrections we neglect. In Ref. [13J, a prescription for re- 
moving such systematics by improving the heavy-quark 
action has been presented. Such improvements are not 
included in this current work and therefore discretiza- 
tion effects may also be partly responsible for the trends 
seen in the tree-level results of Fig. (for the h chan- 
nel, ~20% rise in the central value by cutting the lattice 
spacing in half; ~10% for the 1"). 

IV. CONCLUSIONS 

We have determined the relative contribution of the 
lowest hybrid configuration to the wavefunctions of 
heavy S-wave mesons. Our approach utilizes the non- 
relativistic approximation to lattice QCD, applicable for 
heavy quarkonia where there is a clear separation of 
radial, orbital, and gluonic excitations and the mixing 
among the corresponding configurations is small. 

Tadpole-improved, tree-level results {c B = 1 in Ta- 
bles [^3 and 0) display hybrid configuration admixtures 
at about (0.063) 2 ps 0.4% probability within the T and 
about (0.15) 2 ps 2.3% for J/ip. The corresponding results 
for the pseudoscalar channels are enhanced by « 3 due 
to spin statistics 0, |2j| . 

Although the radiative corrections for the spin- 
dependent term in our heavy-quark Hamiltonian elude 
evaluation, these factors may provide further enhance- 
ment of the configuration mixing. 

Quenching appears to significantly affect the lattice 
spacing results (see Table ITT1 and Ref. [25j), decreasing 
those determined via bottomonium mass splittings (by 
perhaps 10%, see Ref. 9]). This, in turn, affects the 
quark mass extrapolation and the determination of the 
radiative correction, but the direct effect upon the con- 
figuration mixings is not clear. This may only be an- 
swered by repeating the analysis on lattices with dynam- 
ical quarks. 

Inclusion of higher-order terms within the NRQCD 
Hamiltonian has been shown to alter bottomonium hy- 
pcrfinc splittings significantly (by as much as 15% in one 
case |20j). In principle, such terms can cause the same 
configuration mixings as those determined here through 
only the a ■ B term and may be needed to achieve a 
more accurate measure of the hybrid configuration con- 
tent within the T and rj b . The situation for our charmo- 
nium results is more dire. For these systems, the heavy- 
quark expansion can hardly be trusted on our lattices 
(am c < 1) and higher-order corrections may contribute 
as much or more to the mixing. 

These results have implications for theoretical deter- 
minations of certain quarkonium quantities. For ex- 
ample, since the J/ip — * e + e~ decay should not occur 
directly from the hybrid configuration |32|. our results 
suggest about a (0.15) 2 ps 2.3(2)% suppression of the 
T(J/ip — > e + e~) partial width determined via poten- 



tial models; about a (0.063) 2 ps 0.40(4)% suppression for 
T(T — » e + e~). These effects appear to be rather small 
(the corresponding experimental errors lie at about 7% 
and 4% for J/ip and T, respectively [3^) and this anal- 
ysis lends support to their neglect thus far, at least for 
the T, where our lattice heavy-quark expansion is more 
trustworthy. 

The other side to this issue is that of the vector hybrid 
states containing some admixture of the qq vector con- 
figurations and therefore coupling to e + e~. The above 
(tree-level) results would suggest the appearance of a ccg 
hybrid resonance with a partial width of 

T(ccg -> e+e") ps 0.023(2) x T(J/ip -> e+e") 
ps 0.12(1) keV 

somewhere around E cm — M c5g ps 4.4 — 4.7 GeV (the 
lower estimate arises from using the redetermined value 
of the lattice spacing) and a bbg hybrid resonance with 

T(bbg->e + e-) ps 0.0040(4) x T(T e+e - ) 
ps 0.0053(5) keV 

around M bbg ps 10.8 — 11.1 GeV. Any such fcBg-dominated 
resonance should thus be extremely difficult to see in 
e + e _ i?-scans given this very small partial width and 
(most likely) much larger overall width. The ccg- 
dominated state gives a partial width of the same or- 
der of that seen in a resonance within the same mass 
range: V>(4415) with r(e+e") = 0.47(10) keV IHHlEi; 
1/H4430) with r(e+e - ) = 0.390(74) keV from Ref s. Hi 
l37| . However, this hardly qualifies as evidence of a 
hybrid-dominated structure for such a resonance. Our 
charmonium results are plagued with systematic errors 
which are not easily quantified (given the charm quark 
mass and lattice spacings we use, the mixing through 
other terms in the expansion of the heavy-quark Hamil- 
tonian may be even larger). Besides, any such determina- 
tion would not only require a study of the possible decays 
from the hybrid state (determining the total width) , but 
would also need to consider mixings of the ccg config- 
uration with other cc configurations which lie closer in 
mass. The use of lattices with a finer temporal resolu- 
tion may be helpful in resolving multiple states in future 
simulations. 
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